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The light axigluon model of Tavares and Schmaltz is a viable candidate to explain the Tevatron 
tt forward-backward asymmetry. In this paper we present the forward-backward asymmetries for 
66 and cc systems predicted by the light axigluon. Furthermore, we modify this model to include 
flavor changing couplings of axigluons with the SM quarks. We constrain these couplings from the 
available neutral meson mixing data, and investigate their effects on the rare decay — > fj, + fx~ , CP 
violating D — > h + h~ and isospin violating B — s> it"'*' [i + fi~ decays. We show that a light axigluon 
is a possible explanation of the observed CP violation in D — > h + h~ . 

I. INTRODUCTION 



In the Standard Model (SM), the process qq —> tt is symmetric under the exchange of t and t at leading order (LO). 
When next-to- leading order (NLO) processes are included, there is a small forward-backward asymmetry (FBA), of 
i A§. A g — 0.06 ± 0.01 This non-zero and positive asymmetry means that (anti-)top quarks are emitted preferably 

t-H ' in the incoming (anti-)quark direction. In 2011, CDF [j| Q and D0 collaborations Q measured a high FBA in 
tt production from proton-antiproton collisions. The averaged inclusive asymmetry of the two measurements (at 
' (— I ' 5.4/6 -1 ) is Afb = 0.187 ± 0.037 Q. Furthermore, CDF reported a mass dependent asymmetry [3]: 

^ A% B (m t f > 450GeV) = 0.296 ± 0.067 (1) 

<D ■ A% B {m tl < 450GeV) = 0.078 ± 0.054 (2) 

13, ; 

There is another asymmetry in the tt system that is measured at the LHC, the charge asymmetry, that is the average 
' rapidity difference between top and anti-top quark. The SM prediction for this asymmetry is A c = 0.007(1). The 
[ averaged LHC charge asymmetry is Ac = 0.001 ± 0.014, which is close to the SM value Q. 

There are various new physics (NP) models to explain the FB asymmetry, many of which are in tension with the 
0^ ' LHC charge asymmetry, like sign top production, and the tt cross section. In this paper we will consider and modify 
the light axigluon model suggested by Tavares and Schmaltz [3], which is still a viable candidate [13]. This model 
contains axigluons (massive counterparts of gluons) that have a mass below 450 GeV, and that have universal axial 
t-H , couplings to quarks. It gives a considerable tt FBA without significantly changing the tt cross section. It is also 
consistent with the charge asymmetry measured at the LHC Q . 

For details of the model, see [3, EH. To summarize, the model has an extra SU(3) symmetry group, hence the 
gauge symmetry is SU(3)i x SU(3)2 X SU(2)w X U(1)y- Introduced with this extra symmetry group is an extra 
set of up- and down- type quarks, and a scalar field $, which acquires a vacuum expectation value (vev) to break 
SU(3)i x SU(3)2 into the diagonal SU(3) C of the SM. Through this symmetry breaking, one combination of the two 
SU(3) gauge fields acquires a mass. This massive color octet is called the axigluon, and its massless counterpart is the 
SM gluon. Similarly, there are combinations of fermions that become exotic heavy quarks and the SM light quarks, 
allowing the axigluon coupling to the light quarks to be a free parameter. Also, there are no gauge anomalies due to 
cancellations from the additional quarks. Gluons couple to both the SM and heavy quarks with the same strength as 
expected. The lepton sector is exactly the same as the SM, and will not be mentioned throughout this paper. 

In [3], the authors consider only flavor universal couplings of axigluons to the SM quarks. This relies on the strong 
assumption of an underlying global symmetry. This global symmetry is only approximate. Therefore we will not 
assume the existence of an exact global symmetry of the axigluon couplings, which allows flavor changing couplings 
of the axigluons. Other models that have flavor changing color-octet couplings have been proposed in the literature 
[l2T ITjjj. A general Lagrangian with flavor violating axigluon interactions contains the following terms: 



C D Ui^A^gijij + (e2)y)P L ttj + d^A^g^ + (ei) ii )P i d j 

+ ua^{g l uR 8 l3 + {e u R ) lj )P R u 3 + da„A"(g l dR 6 )P R dj (3) 



Here A^ is the axigluon, Ui and di are SM up- and down- type quarks respectively (i = 1, 2, 3 is the generation index), 
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and g % are flavor independent couplings. Color and spinor indices are suppressed for simplicity. The complex matrices 

( n L,R L,R \ / n L R L.R \ 

L,R* n L.R u I L R* n L > R I I A\ 

9ds 9bs > e L,fl = 5„c 5ct (4) 

L.R* L,R* r, I \ L.R* L.R* n / 

9 db 9bs / \9ut 9ct / 

contain off diagonal axigluon couplings of up- and down-quarks respectively. This mixing follows from flavor symmetry 
breaking of heavy and light quarks. These flavor changing neutral currents (FCNC), which can occur at tree level, 
can have interesting effects on FCNC observables. 

The rest of the paper is organized as follows: In Section 1, we investigate the contribution of light axigluons to bb 
and cc FBAs. In Section 2, we constrain the flavor changing axigluon couplings from neutral meson mixing data. In 
Section 3, we investigate the contribution of the constrained flavor changing axigluon model to the following decays: 
B° s -> fjt+fi-, D° -> h+h~, and B -> K^fi+fT. 



II. FORWARD-BACKWARD ASYMMETRIES 

In the Tevatron experiments CDF and D0 , ti production happens mostly via qq — > ti (Fig [T]). 
amplitude of the process qq — > ti including the axigluon contribution is calculated as follows [3, LLL| : 



The square 



\M\ 2 = Ngt(l 
-Ngt 

+ Ng 4 s 



c 2 + 4m 2 ) 
2s(-s + Ml) 



(-s + M\Y + 



Y\M\ 



P{-s + M\f 



{{-s + Mlf+T\Mlf 



[9 q v 2 



-c 2 +4m 2 ) + 2g t A g q A c 
gf] + c 2 + 4m 2 ) + g t2 {l + c 2 - 4m 2 )] + 8gv9A^9 9 A c 



(5) 



where N — | is the color sum, s = — (pi + P2) 2 is the partonic total momentum, /3 = y 1 — -j 1 = y/T— 4m? is the 
velocity of the top quark, c = j3 cos 9 where 6 is the angle between the incoming q- and the outgoing t-quark, Ma is 
the axigluon mass, and Ta is the width. Vector and axial couplings of the axigluon are defined as: 



9v 



Q 1 Q 
q _ 9R+91 



Q 9r 

9a = — 



9l 



In Eq. [5] the first term comes from the SM gluon exchange, the second term is the interference between the gluon 
and the axigluon channels, and the third term is the axigluon s-channel (See FigJT]). The FBA comes from the terms 
that are proportional to the odd powers of cos 8, since cos 9 is odd under 9 — > tt — 9. In order to accommodate the 
measured ti FBA, the axigluon should give a large FBA without affecting the ti cross-section, which is close to its 
SM value. A light axigluon (Ma = 100 — 400 GeV) with a large width Ta — O.IMa, is shown to agree with both the 
ti FBA and the ti cross section [H, [l(| • Assuming approximate parity symmetry, g l v needs to be small. We will, as in 
@, take g l v =0. 




g^(g q B PR + glPL) 



ma 




FIG. 1: Diagrams contributing to qq — > tt amplitude. Left diagram is from the SM gluon exchange, right diagram is 

the axigluon exchange. 
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We will define the FBA, Afb, through the forward and backward scattering cross-sections: 



a+ = cr(0 < < tt/2) 
cr_ = a(w/2 < 9 < 7r) 

where a — J dcos9 \ A4\ 2 is the total cross-section. Then the FBA is 



(7 + — (7 



(6) 



Rn 2 ~ S + M A.. 

PyA(-s+MlY+T\M\ 



(7) 



5j(l + £ + 4 ™ 2 ) + 



s(-i+Mj|)2 



T3l(l + f -4m 2 ) 



2 [(_s+Ml) 3 +r2Mi] 



where <^ = <7^ = qa ■ This FBA is plotted in Fig[3] as a function of the tt invariant mass m t i . The values for coupling 
constants for each axigluon mass (100, 200, 300, 400 GeV) are taken from |l0|, where the authors performed fits to 
both the Tevatron FBA and the LHC charge asymmetry, and chose the coupling constants that would best fit both 
of them. 

An obvious prediction of the light axigluon model is FBAs for bottom and charm quarks. The graphs for these 
asymmetries can be seen in Fig|4"aland Fig|4b[ Since both m& = 4.2 GeV and m c — 1.27 GeV are very small compared 
to the predicted axigluon mass (100s GeVs), the FBA structure is almost the same for both quarks. As can be seen 
from FigHl the predicted FBAs are quite large even for low energies, m b g ~ 200 GeV, so the search needs not to go 
to high invariant masses. Furthermore, the crossing at m q q — Ma (<Z = b, c) is expected to be clearer compared to tt 
production, since Ma >> 2m q . Consequently, measuring the bb FBA will be a better way to find out the axigluon 
mass. Measuring the bb FBA is also suggested in 

III. CONSTRAINTS ON FLAVOR CHANGING AXIGLUON CURRENTS FROM NEUTRAL MESON 

MIXINGS 

In this section we will use the data available from neutral meson mixings (B® — B^,K° — K°,D° — D° mixing) 
to constrain the FC axigluon coupling matrices, e£ R and t d L R , in Eq@] Neutral mesons "mix" because the flavor 
eigenstates (M°, M°) of the SM Hamiltonian are not the actual mass eigenstates (Ml^h, L and H for light and 
heavy respectively). The mass difference between Ml and Mh will be one of the constraints we will use for each 
mixing. Since these processes are FCNC processes they happen via loop diagrams in the SM, like electroweak (EW) 
box diagrams in Fig{6j These diagrams are suppressed by Cabibbo-Kobayashi-Maskawa (CKM) matrix elements due 
to flavor changing EW interactions. However, with flavor changing axigluon couplings we have neutral meson mixing 
at tree level. Therefore the FC couplings are constrained by the mixing data for B® — B® and K° — K° and D° — D a 
mixing. 

For the following calculations we take Ma — 400 GeV, which gives us the least stringent constraints. We explain 
the methods we use to constrain each coupling constant in the respective sections. In each section we try to give 
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FIG. 2: tt forward backward asymmetry, Apg vs tt invariant mass m t -. Afb only changes sign when Ma > 2m t . 

The values for the coupling constant gA are taken from [Toj . 
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(b) CDF graph of Ap g vs tt invariant mass m t j [f| . 



FIG. 3: Comparison of the light axigluon prediction and the recent CDF measurement of tt Afb- Due to the large 
error in the data, Ma=400 GeV is not precluded. The values for the coupling constant <?a are taken from 10]. 
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FIG. 4: Forward-backward asymmetries for bottom and charm quarks vs qq invariant mass m q q, for q = b,c. All 
asymmetries cross zero when m q q = Ma which is much higher than m q q = 2m q . The values for the coupling 

constant gA are taken from [lOj. 



the most general forms that can be used to constrain these couplings, but at the end we assume axial couplings for 
simplicity. We prefer to give constraints on real and imaginary parts of gfj = quark flavors) rather than gij itself, 
since the meson mixing amplitudes that we use for the constraints involve the square of the coupling constants. One 
can find the constraints on the real and imaginary parts of the coupling constants themselves, assuming neither of 
them are zero. The magnitudes of the coupling constants are constrained using \z\ 2 = |z 2 |, which is true for any 
complex number z. A summary of the results can be seen in Table. [J. 



Coupling 


Constraint 




\9ij\ 


Qbd 


5° - B° 


< 2.10 x 10" 7 


< 4.58 x 10" 4 


gbs 


&s — ts a 


< 2.55 x 1CT 6 


< 1.83 x 1(T 3 


gds 


K° - K° 


< 6.13 x 1(T 13 


< 3.11 x 1CT 5 


guc 


D°-D° 


< 4.89 x 10" 9 


< 1.47 x 1(T 4 



TABLE I: Origin of constraints and upper bounds on the imaginary part and the magnitude of FC axigluon 

couplings 
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FIG. 5: Integrated asymmetries from Fig. [4] 
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FIG. 6: Electroweak box diagrams that contribute to B a — B a mixing 



A. B° q - Bq Mixing 



The SM and the axigluon contribution to — mixing amplitude, where q = s,d for B° s and B° respectively, 
can be seen in FigjTJ The gray box in the first diagram is the sum of all EW box diagrams that have u, c, ^-quarks in 
the loop. However, for B® — B® mixing, the Hoop is the most important one. The SM contribution, including NLO 
corrections, to this amplitude is given in [l"> as: 



M 



SM 



W 



16tt 2 



a s (pb) 



4tt 



(bq)v~A(bq)v-A 



(8) 



where Gp is the Fermi constant, My/ — 80.4 GeV is the W-boson mass, are the CKM matrix elements that mix 
i— and j— type quarks, a s (/if,) is the strong structure constant evaluated at /if, ~ 0(m,b), and J/ (~ 1.627 for / = 5, 
/ being the number of active flavors at the mixing scale) is a constant that comes from the running of the coupling 



r{gtPL + gtP R ) 



FIG. 7: The SM and the axigluon contribution to B q — B q mixing 
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coefficients. Also (bq)v-A = &7 M (1 — 75)9, and 

So(x t ) 



4x t - llx? + x? 



3x? 



4(1 ~x t f 



2(1 -xtY 



In x t 



(9) 



with Xj = -^3-- This function iSo(x) is one of many similar functions called Inami-Lin functions. They are loop 

functions that arise in box and penguin diagrams in the SM and are calculated by Inami and Lin in 16]. For 
m t = 172.1 GeV and M w = 80.4 GeV, S (x t ) ^ 2.51. Thus the SM mixing amplitude is: 



M 



SM 



' 16tt 2 



(2.U)(V t * b V tq ) 2 (bq) V - A (bq) V - A 



(10) 



The second diagram in Fig 13 is the axigluon contribution to the mixing amplitude. This part can be written as: 



M n 



2C"(/i) < g& (bi^PLTfjqj) (b kl »P R TZ iqi ) 



(11) 



where C(/z) and C'(/z) are coefficients that run with the scale /«, and Pl,r = ^(1 T75) are projection operators. The 
first two terms in the square-brackets are similar to SM (V — A) ® (V — A) couplings, but the last one is of type 
(V — A) <E> (V + A) which mixes left- and right-handed quarks. By using the following sum rule, 



1 ij 1 kl — 



and gamma matrix identities, 



we can write Eq[TT]as follows: 



(i^PLh(rp L )ki = -(^plU^pl)^ 

(rfMui-fPRiu = 2(p R ) a (PL) kj 



(12) 



(13) 
(14) 



Max = 



Ci 



{byi'PLq)(h' 1 PLq) 



9bq9bq 



(by»P L q)(lr/i i P R q) + C 2 {bP L q)(bP R q) 



(15) 



where N c is the number of colors, and Co, C\ and C 2 are the renormalization group (RG) evolved coefficients for the 
corresponding 4-quark operators. We take Co from the SM since those operators also exist in the SM. On the other 
hand, we calculate the RG evolution of C\ and C2 here. Before doing that, let us define the following quantity: 



(B° Q \n 



SM 



m\u SM \m 



(16) 



in order to compare the new physics (NP) and the SM contributions to the mixing process. We need the following 
matrix elements 11211: 



{B° q \b^{\ ± 75 W(l ± J 5 )q\B o ) = ~m BO flB g 



(B° q \br(l + 75)<Z&7 M (1 - 75)<Z|5 ? °) = - 
(B° q \b a r(l + 75)9^/37^(1 " lMB° q ) 



™ a J B B RL 

i 

7 

3 



m B olB q B 2q 



where m e0 is the B q meson mass, fs q is the decay constant, and B the bag parameter (B q 
B-meson decays). These matrix elements come from lattice calculations. Recent lattice-calculated values for these 



B lq 



(17a) 
(17b) 
(17c) 



B*? ~ 1 for 
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constants can be found in [T3, Gil and references therein. Using Eq[T7] in EqflO] and [15] we can write the SM and 
axigluon matrix elements for the mixing as follows: 



r<2 yr2 

(B° q \H SM \B q } =^-f(2.lA){V; b V tq fm B jl q (18) 



1 



(B v \U ax \B ql , 



M A 



^ (CO 2 + c<) 2 ) + (lr - 7 -§) « 



m HO fl (19) 



Now let us look more into the coefficients Ci(/i) and C2(//). They obey the following RG flow equation [191 ] : 

<?(m) - fO^O*) (20) 

where 
and 

is the anomalous dimension matrix. Here 7 is the LO and 7 1 is the NLO contribution to the anomalous dimension: 



/ 6 19 \ / ML A- .15 _ 22_ f 200 at 6 44 f 

_ I JVe iZ I 1 _ I 6 T 3 c N c 3 J 

I -tN +j-J' 2 I n 7Vc + j__ 2 / -"SjV* + 479 + « + *tf c /-JL/ 



Eql2"01 defines how to run Ci and C2 down from a higher scale (Ma ~ 400 GeV) to a lower scale (/1 ~ 5 GeV) . In 
this paper we only take the LO evolution of the axigluon coefficients. We can diagonalize 7 in Eqj20]to get 

_ 6/JV c _ 6/N c -6N c 

a(n) \ ^0 2 / a(/i) \ ^0 



c ^ = K(^m>) -k{*m7>) (21) 

_ 6/W c -6iV c 



C2 ^ - I ^7) ) ( 22 ^ 



where j3 = 11N f 2f . Also, the SM-like coefficient C runs as follow pj|: 



W=(-^V 2 " (23) 



v a(M A ) 

For — B° mixing, N c = 3, / = 5, and ^ ~ 5 GeV. With these, the axigluon matrix element becomes: 

m o/s 

(B°|^ax|B°) * (-0.18(( 5hg ) 2 + ( ff «) 2 ) - 0.73 gfa*) (24) 



Thus EqfT6] reads: 

12rr 2 ( -0.18((g^) 2 + «) 2 ) - 0.7 3 gfcgg) 
(%Af^M»(2.14)(V£V t ,) 

As in the flavor conserving part of the light axigluon model, we assume axial couplings: g^ q = —g^ q = fjbq- 



A =1 + 12Tr 2 (0.37)( g&g ) 2 

(2.14)0^^1(^)2 ^ 

In [2I] one can find Re(A g ) and Im(A 9 ) for B® and B° mixing. In the next two subsections we are going to look at 
both cases separately. 
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1. B° — mixing 



In this paper we use the CKM basis given in [21(. In this basis, the relevant CKM matrix elements, Vtb and Vt s , are 
real, which makes the SM contribution real. Hence, the real and imaginary parts of g% g can be constrained separately 
from the real and imaginary parts of A s . From [20[, we take the following boundaries: 

0.85 < Re(A s ) < 1.27, |Im(A s )| < 0.32 (27) 



To be conservative, all parameters are taken at the 3cr boundaries of the fits from [20], since there is some tension 
with the SM at 2a. Using these values and Eqf26]with q = s, we get the following constraints: 

|Re( ff 2 s )| < 2.15 x 10- 6 , \Mgl)\ < 2.55 x HT 6 (28a) 
=► \g bs \ < 1.83 x 10- 3 (28b) 

2. B° - B° mixing 

For B° — B° mixing, one of the relevant CKM matrix elements, Vtd, is complex, therefore the real and imaginary 
parts of g"l d can't be constrained separately. Let us work this through in more detail. We can write Eqf2S]as follows: 

A d = 1 + Cg 2 bd (29) 

where 

12tt 2 (0.37) 

6 - {2.u)GiM* w Mi{v: b v td y lM + L72i 

Notice that the real and imaginary parts of C are very similar, so we will take Re(C) ~ Im(C) = a. Then the 
constraint equations are 

|Re(A d ) - 1| = a|(Re( 5 2 d ) - Im( 5 2 d ))| (30a) 

|Im(A d )| = a|(Re( 9 2 d ) + Im^ 2 ,))! (30b) 

Again from (2^, we take the following bounds (at 3a): 

0.62 < Re(A d ) < 1.36, -0.39 < Im(A d ) < -0.02 (31) 

Note that, for the Im(Ad), the SM value of is slightly outside the 3tr allowed region, but we choose to disregard this 
discrepancy in setting the limits, since it is very small. Now, notice that |Re(A,j) — 1| < 0.38, and |Im(Ad)| < 0.39. 
Putting these all together with Eg I3TJI we get similar constraints for Re(<? 2 d ) and Im(g 2 d ) such that g\ d lies in an 
approximate circle of radius ^0.39/a. So we get 

|Re( 5 2 d )| < 2.10 x 10- 7 , |Im(.g 2 rf )| < 2.10 x lO" 7 (32a) 
=> \g bd \ < 4.58 x 10- 4 (32b) 

B. K° - K° Mixing 

The last constraint on e d comes from K° — K° mixing. The SM and the axigluon contributions to K° — K° mixing 
are shown in Fig. [8] The LO SM contribution again comes from EW box diagrams, but in this case c- and i-quark 
loops both need to be considered. The NLO mixing amplitude from the SM is [jfjj : 



.G%M? 



M SM = -i f^f (X 2 c (lM)S (x c ) + X 2 t (Q-57A)So(x t ) + 2X c Xt(0A7)So(x c ,x t )) (a s (v))- 2/9 



1 + ^T J3 



(sd)v~A(sd) v ^ A 
(33) 

where 



x t Sxt 3x 2 

~x c ~ 4(1 -a; t ) ~ 4(1 -a;*) 2 



S {x c ,xt) = Xc In — r - -j- —\nxt) (34) 
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FIG. 8: The SM and the axigluon contribution to K — K mixing 



is another Inami-Lin function, A; = V* s Vid, J3 — 1.895, and /j, ~ 0(1 GeV). The relevant matrix elements for kaon 
mixing can be found in [22]: 

(K°\sr(l± l5 )ds^(l± l5 )d\K°) = \m K J 2 K B K (35a) 
(^ |^(l+ 75 )^(l-75)rf|^ > = (7.8)^/1^ (35b) 
(K°\s a ^(l + 75)^^7^(1 - j 5 )d a \K°) = (30A)^m K J 2 K B K (35c) 
Using EqESJin Eq|33]we get 

c 2 m 2 

(K°\-H SM \K ) = ^-fm K J 2 K B K [A*(0.001) + Af(1.776) + 2A c A t (0.001)] (36) 

The axigluon contribution to kaon mixing can be written using Eq |15l with the substitution of appropriate quark 
operators. Then using Eg 1351 the axigluon matrix element becomes: 

(K°\n ax \K°) = m *°^f K {~^((9 L ds ) 2 + (9f s f) + (-1-73C! + 10.13C 2 )^«} (37) 

The RG evolution of C\ and C2 is similar to the one in the previous section, with a difference in scales. This time we 
need to run down to \i ~ 0(1 GeV) (and hence the number of active flavors will be / = 3). After running down the 
coefficients Co, C\ and C2 to /i ~ 0(1 GeV), the axigluon matrix element becomes: 

(K°\H ax \K°) = mK "ljf K {-0-16((^J 2 + ( ff £) 2 ) + 48.68 g L ds g*} (38) 
Assuming again axial couplings: g% s = — g ds = gds, the axigluon matrix element is: 

(K°\n ax \K) = m «^ J? * (48.36)gg, (39) 

Now, following a common notation [21], we define 

M 12 = (K°\H eff \K ) (40) 



This matrix element M\% is the off-diagonal element of the "mass matrix" in the full Hamiltonian H — M + |T, and 
W e ff is the effective Hamiltonian, that includes both the SM and the NP interactions, for the mixing process. The 
off-diagonal elements of M are related to the mass difference between heavy and light mesons, and the off-diagonal 
elements of T are related to the decay of these mesons. The interested reader should refer to 21] and references 
therein for more information on meson mixings and decays. The CP violation in meson mixings comes from a possible 
phase difference between Myi and Ti2, which depends only on the short distance part of the matrix element Myi. 
The long distance interactions, which come from on-shell states in the loops (Fig|3]), are CP conserving. The axigluon 
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does not contribute to the long distance part of this amplitude at LO. The mass difference, Am, between heavy and 
light mesons is 

Am = 2\M 12 \ (41) 

Notice that Am gets affected by both the short distance and the long distance parts of the effective Hamiltonian. 
Unfortunately, calculation of the long distance contributions is difficult [15| . In this paper we assume that long 
distance contributions are at most 50% of the total mass difference, hence My, M = 2M§f. 

In order to constrain the coupling constant gd s , we use the mass difference [2lj | and the imaginary part of Mi 2 22}: 

Im(M^ p ) = (1.7 ± 1.6) x 10~ 18 GeV (42) 
Am = (3.483 ± 0.006) x 10~ 15 GeV (43) 

Consequently, we get the following constraints: 

\Re(gj s )\ < 9.64 x 10~ 10 , |Im(«4)| < 6.13 x 10" 13 (44a) 
=>■ \g ds \ < 3.11 x 10- 5 (44b) 

Finally, we can write e d constraints as follows: 

9.88xl0~ 9 2.08xl0~ 4 \ / 0.311 4.58 

|Im(e d )| < I 9.88xl0~ 9 7.69xl0- 4 , \e d \ < 0.311 18.3 | x 10~ 4 (45) 

2.08xl0~ 4 7.69xl0~ 4 / \ 4.58 18.3 

where Im(e d ) is calculated assuming Im(gij) < Ke(gij) and Im(<?jj), Ke(gij) 7^ 0. 



C. D° - D° Mixing 



D° — D° mixing, like the other meson mixings, happens via EW box diagrams in the SM, shown in FigJHl Calculation 
of the LO SM contribution to the mixing amplitude is similar to K° — K° mixing where c— and t— loops are changed 
with s- and 6-loops. D° — D° mixing suffers more from long distance contributions as compared to K° and B° mixing. 
Consequently, NLO corrections to the mixing amplitude are not calculated in the literature. In this section, we only 
consider the LO short distance amplitude, and assume that long distance effects are at most the same order as short 
distance ones. The LO short distance SM contribution to the mixing amplitude (at /i = 2 GeV) is: 



C 2 M 2 

M SM = -i ^^ (0.80) {X 2 s So(x s ) + X 2 b S (x b ) + 2X s \ b S (x s ,x b )} (uc) v . A (uc) v . A 
C 2 M 2 

- ~i h 7 i x2 s ( L24 x 10_6 ) + X l ( 2 - 18 x 10_3 ) + 2X - X b ( 9 - 23 x 10_6 )} {uc) V -a(uc)v-a 

lD7T z 



(46) 
(47) 



u 



+ 



it 



u 



FIG. 9: The SM and the axigluon contribution to D° — D° mixing 
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For the matrix elements, we assume similar relations to Eq[l7] since m D0 ~ m c + m u ~ 1.86 GeV. The axigluon 
contribution can be written using Eq |15l with the substitution of the appropriate quark operators: 

(D a \H S M\D Q ) ~= ^ w {A 2 (1.24 x lO" 6 ) + A 2 (2.18 x lO" 3 ) + 2A s A b (9.23 x 10~ 6 )} m D of DO (48) 
(D°\H ax \D°) |-0.18(( 5 £) 2 + (.g,i) 2 ) - lMf&fa} m D of DO (49) 



Once again, we assume axial couplings g^ c = ~9uc = 9u 



[irj{ . Z))i . = !^ (0 . 68)£c (50) 



Constraints on g uc come from the D° — D° mass difference, Amu = 1.57 x 10 14 GeV [2l[, and the ratio v/p, where 
q and p are the coefficients that describe the flavor eigenstates D°, D° in terms of mass eigenstates D° Hl D Q L : 

\D H ) =p\D°)-q\D°) 
\Dh) — p\D°) + q\D°) 

One can show 



q = M* 2 - fr* 2 



In L>° - D° system, T 12 ~ M12 [H|, and so we have 



q = 2M^ 



(52) 



We use Eq |52l to constrain the imaginary part of g uc . The real part is already more constrained from the mass 
difference. From [2lJ, we take the following values for the magnitude and the argument of ijp\ 

= 0.60 Arg ( - ) = -22.1° (53) 



P 

The constraints on g uc from these values are as follows: 



\Re(gl c )\ < 2.09 x 10" 8 , |Im(. 9 2 c )| < 4.89 x 10~ 9 (54a) 
=> \g uc \ < 1-47 x 10- 4 (54b) 

Unfortunately, there are no other mesons with which we can investigate the up-sector further. However, the neutral 
D-meson system has other interesting features, like the CP asymmetry in D° — s- h + h~ decays that was measured in 
2011 [23j]. We look more into the contribution of FC axigluons to this process in the next section. 



IV. SOME EXAMPLES OF FLAVOR CHANGING AXIGLUON CONTRIBUTIONS TO MESON 

DECAYS 



In this chapter, we will check the effects of the FC couplings on several SM processes, namely the rare decay 
B° s — > the CP asymmetry in D° — s- h + h~ decays, and the isospin violation in B ^ — > K^fi + fi^ decays. We 

do not expect a significant contribution to B® — > decay from the axigluons, since it is affected through axigluon- 

penguin loops at LO. However, when FC axigluon currents exist, processes like D° — > h + h~ and B ^ — > K^*'^, + fj," 
can happen at tree level. Therefore one would expect to get an appreciable contribution from axigluon induced 
channels. These decays are chosen because they are of current experimental interest [23T426T ] . 
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A. Rare decay B a 



In the SM, this decay is predicted to happen very rarely, with a branching ratio of (3.5 ± 0.30) x 10 -9 [27] . This 
is very close the recently measured branching ratio of 3.2+^2 x 10 -9 by LHCb [24| HBj]. These new results constrain 
the NP one can have that would affect this branching ratio. As we will see in this chapter, axigluon contribution to 
this decay amplitude is at least 3 orders of magnitude smaller than the SM contribution. Consequently, the measured 
branching ratio data does not impose further constraints on the coupling constant g bs . 

This branching ratio is so small in the SM as it occurs through EW penguin and box diagrams. The SM effective 
Hamiltonian, including NLO corrections, is calculated in [15] as: 



where 



Y(x f ) 
Y (x) 

Yi(x) 



Y (xt) - 
'4- 



A* A* 



.G, 



a 



y/2 2tt sin 2 6 W 



-{V: b V ts )Y(x t ){sb) v „ A {w)v-A 



4?r 



4a; 



1 - x 
- 16a; 2 



Yx{x t ) 

3.c 
(1^)- 
ix 3 Ax 



In a; 



10a; 2 



3(1 -x) 2 
2x + x 3 
' (1 - x) 2 



L 2 {1 -x) + 8x 



(1-x) 3 
dY (x) 



■ In a; 



2x - 14x 2 



2(1 - x) s 



dx 



In a;,. 



• In 2 x 



(55) 

(56) 
(57) 

(58) 



and 0w is the Weinberg angle. In Eq|58l x^ = where /u, is the renormalization scale (/x 
GeV and M w = 80.4 GeV, Y(x t ) ^ 0(1), and EqlSUcan be written as: 



Hsm{B s 



a 



A* A* 



y/2 2tt sin 2 9 



0(m B )). Form t = 172.1 



(59) 




A* 



FIG. 10: The axigluon contribution to — > ji + ji 



The axigluon contribution to this decay is shown in FigfTUl which is a penguin loop with an axigluon. This is the 
LO contribution, since the axigluon does not couple to leptons. In this diagram, b— and s— loops are more important 
than the d— loop, since the d— loop is doubly suppressed by the small coupling constant. We only give the 6-loop 
calculation, the calculation of the s-loop follows simply by substituting an s-quark wherever there is a 6-quark. For 
simplicity, we set the external momentum to zero. The amplitude is 



Max = 



9Us 



4 cos 2 OwM"; 



-(A*At) 



V-A 



d n l 



s^N^j^b 



(27r)« (l 2 + m 2 ) 2 (l 2 



M 2 A ) 



and 



N u = (g^P R + gb s P L )(-t + m 6 )7>& - a bl5 )(-t + m b ) l5 



{9bsPR + 9bsPL)h v t{v b l5 + a bl5 ) + m 2 {ggP R + g^P L ) j > fc75 - a b ) 



(60) 



(61) 



where Vf = T 3 — 2Q f sin 9w and a,f — , T3 and Q j being the isospin and the electric charge of the fermion / 
respectively. We do not keep the terms proportional to m b , since they integrate to zero under the loop integration. 
Using some gamma matrix identities we can write (in n = 4 dimensions) 

rN» lfl = -fb&PR + gt.PL) [l 2 (-VbK + a b ) - 2m 2 (v bl5 + a b )] (62) 
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There are two integrals to deal with. One of them is the following divergent integral: 

d n l I 2 



h = 



(27r)» (Z 2 +m 2 ) 2 (Z 2 + M 2 ) 



d n l I 2 
2 / dXldX2 I (27r)« [I 2 + ( Xl + x 2 )m 2 + (1 - Xl - x 2 )M\f (63) 



Using dimensional regularization for n = 4 — e, we get 



h (L) 2 J dXldX2T Q {(x 1 +x 2 )m 2 + (l-x 1 -x 2 )Ml) 



f d Xl dx 2 (\n (- ^ —f) + const +- + 0(e) ) i (S4 ) 



(4tt) 2 y Z V U^i + ^W + ti-n-^M] 



r 



This divergence arises because the part of the theory we are considering is not complete, for example there are 
additional (heavy) quarks in the full theory, and their inclusion should lead to cancellations of these divergences. The 
full theory is renormalizable [9j , so we do not worry about the terms that can be canceled through the short distance 
contributions of the heavy quarks. 



h = , . / dx\ / dx 2 In 



(4tt) 2 Jo Jo V Oi + x 2 )m 2 + (1 - xi - x 2 )M 2 A 

ln T7T ~ 77^2 / dxi / dx2 M 1 _ x i - x 2 + ( x i + x 2 )y b ] (65) 



(4tt) 2 \M\ J (4tt 



Bi{y b 

where y b = jjjfr and 



y(y-2 ) 

4 ' 2(y-l) ' 2(y-l) 



J? i(^) = -7 + TT7-^T+ r M 2 ln ^ ( 66 ) 



The first term in EqJSS]can be used as an estimate of the short distance behavior when fj, ~ A ~ O(Ma), where A is 
some cut-off scale for NP. Consequently, we take the integral to be 



(47T) 

The other integral is finite: 

f d n l 1 f f d n l 



^-8^1 (I/) (67) 



(2tt)« (Z 2 + m 2 ) 2 (Z 2 +Af 2 ) 7 1 ' J (2tt)« [J 2 + (an + x 2 )m 2 + (1 - Xl - x 2 )M 2 A Y 

^ dx\dx 2 - r — (68) 



(4tt) 2 M 2 J " 1 - Xl - x 2 + { Xl + x 2 )y b 



B 2 (Vb) 



where 



1 lny 
B 2 {y) = - —2 

Inserting Eq|67| and [68] into EqipTJl we get 

M ax = -^^(w)v-A s [-f(g&P R + 9tsPL)[(-v b l5 + a b )Bi(y b ) + y b (v bl5 + a b )B 2 (y b )}] b (69) 
Next, realizing that yB 2 (y) << B±(y) for y << 1, we neglect the part ^B 2 (y), and write as an 0(1) estimate 

M a , ^ -^-{B x {y b ) + Bi(y s )){nfi) v ^ A [s^(V - A l5 )b] (70) 
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where 



V _ gj\+9bs 



and 



fJl. 



together with 
B® — > decay hamiltonian becomes 



V = gY s a b - gtvb 
A = gl s vb - gt a b 

In this paper, we take g^ s = 0. Thus the axigluon contribution to the 



Uax = i — ( B i(Vb) + Bi(y a ))(sb) V -A [w)v-A 



Now we can compare Eq[7T]and EqUH for Ma = 400 GeV (at M w scale): 



(H 



(H 



0.001 



SM/ 



(71) 



(72) 



The ratio gets smaller for smaller axigluon masses. Since the axigluon loop contributions are very small compared to 
the SM, the uncertainties in the calculations should not be a worry. Hence flavor changing axigluon couplings under 
already considered constraints do not affect the B® — > branching ratio in a noticeable way, and so this decay 

does not give further constraints on the coupling constant g os . 

B. CP Violation in D° -¥ h + h~ decays 




U. 




u 



FIG. 11: CP violation in D° 



h+h- 



The CP asymmetry in D° —> h + h decays is defined as follows: 

T(D"^h+h-)-T(D"^h+h-) 
CP T(D° -> h+h-)+T(D° -> h+h-) [ ' 

In 2011, the LHCb measured this asymmetry to be ~1%. There still appears to be much confusion about the origin 
of this asymmetry in the SM (28l - [30l | . The long distance effects in D-decays make the calculation of relevant hadronic 
matrix elements very difficult. Furthermore, the charm quark might not be heavy enough to trust perturbation theory 
at this scale. Naively, one would expect the CP violation at LO to come from decay channels like the first diagram in 
FigHU which is estimated in [H 111 to S ive ~0.1%. The discrepancy between 1% (experiment) and 0.1% (SM) might 
be due to new physics. However it also might be contained in the SM if some matrix elements of penguin operators 
are much larger than the estimates given by dimensional analysis (28j . In this section, we compare the contributions 
of the two diagrams in Fig[TT] to the CP asymmetry. 

The effective Hamiltonian that comes from the SM diagram in Fig 111] can be written as follows [311 ] : 

G 



where 



G 

Hsm = {C 3/5 (uc) V - A (d^d) + C 4/6 (u Q c ;3 )y_ yl (d ( 37^d Q )} (74) 



C 3/5 (M W ) = -lc 4/6 (M w ) = - a '^ Eo(x b ) (75) 
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and 



E (x) = -- H x) + { m „ x)3 + 6(1 „, )4 M*) " 3 (76) 
is another Inami-Lin function. The subscripts of the coefficients and the choice of writing the quark operators in this 



way is a slight variation of what Buras does in his paper [31( . He gathers operators and their coefficients in a way 



that is easier to keep track of in the RG flow equations. Here we do not RG flow the coefficients, and compare the 
SM and the axigluon parts at My/ ~ Ma — 0(100 GeV). The effective axigluon Hamiltonian that can be written 
from the diagram on the right in Fig jlll is: 

9uc9s f l/_ w w j ^ , 1, 

where q s ~ % @. Comparing EqEZ] with EqEl for M A = 400 GeV, we see that the upper bound for axigluon 
contribution is an order of magnitude larger than the SM contribution. This upper bound grows with decreasing 
axigluon mass, to ~40 times the SM contribution for Ma = 100 GeV. This result can explain the measured CP 
asymmetry at the LHCb without requiring the matrix elements to be larger than given by dimensional estimates [23| . 



C. Isospin violation inB^ K^*' fi + \i decays 



In the SM, B — > Kfj, + n~ decay follows, at LO, from EW penguin and box diagrams that do not involve the d(u)- 
quark, which is then called the spectator quark (Fig ll2[) . For this decay, an observable, the isospin asymmetry Ai, 
can be defined as follows: 



Ai 



) 



T(B° -> K Q n+ii-) + T(B+ -> K+fi+fi-) 



(78) 



where T(i /) is the partial width of the i —¥ f decay. A similar asymmetry is defined also for B — > K* . 
We can see that in the spectator quark approximation, this asymmetry is zero, since there is no difference between 
the decay of the neutral and charged B— meson. If we consider diagrams in which the final pair is emitted 

from the spectator quark (FigfT3|). there would be a non-zero isospin asymmetry due to the different charges of the 
spectator quarks involved in neutral and charged B— meson decays. In the SM, the asymmetry for B — > K*[i + [i~ is 
expected to be around —1% (32l.l33j. Although there is no clear prediction for the isospin asymmetry in B —> Kfi + 
from the SM, one might expect it to be similarly small, almost zero (26[. The isospin asymmetry that is measured 
at the LHCb is consistent with the SM for the B — > K*/j, + fi~, however it deviates from zero with 4.4ct significance 
for B — s- K/i + fj,~ [26| . The SM prediction might be enhanced by more precise hadronic matrix element calculations, 
however there might as well be NP involved in these decays, like flavor changing axig luons (Fig|l3b|. 





u, c, t 



d,u d.u 




u, c, t 



d, u d,\ 



FIG. 12: Isospin conserving diagrams in the SM for B 01 -^ —> K(*>/j, + ii decays 



In order to compare the axigluon contribution to this isospin violating process with the SM contribution, we assume 
that the EW penguin diagram in Fig ll3al is as important as any other isospin violating diagram in the SM, if not the 
most important one. Therefore, instead of performing comprehensive calculations of the SM contributions, we only 
compare the two diagrams that are shown in Fig ll3l Furthermore, we only consider the parts of these diagrams that 
are responsible for B — > K° decay, since the emission of the final muons are the same in both cases. 

The SM contribution to the hadronic part of the effective amplitude from Fig |13al can be written as follows: 

M SM - X~^E (x t )(bs)v~A(dd) (79) 
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w 




(a) One of the isospin violating SM diagrams (b) Axigluon contribution to isospin violation 



FIG. 13: Isospin violation in _B°( + ) — > K°^fi + fi decay 



where i — u, c, t. 

The axigluon part of the same amplitude from Fig |13bl is: 

9bs9s 



M ax * ^(bs)v-A(dd) (80) 



6Mj 



Comparing EqEHland EqlSDlfor M A = 400 GeV (at M w scale), we get 



(Ho 



0.3 



(Hsm) 

The axigluon contribution is at most the same order as the SM one when Ma = 100 GeV. 



(81) 



V. CONCLUSION 



The light axigluon model of Tavares and Schmaltz is an experimentally allowed modification of the SM and a viable 
explanation of the CDF ti forward-backward asymmetry P, [13] ■ In this paper we used this axigluon model to predict 
bb and cc forward-backward asymmetries. They are expected to be large and depend on the invariant mass of the 
quark pair. This mass dependence is a useful tool to investigate the mass of the axigluon. 

We also modified this flavor conserving axigluon model to include flavor violating couplings between the axigluon 
and the SM quarks. These couplings are constrained by neutral meson mixings, and the upper bounds on their 
magnitudes are in the range 10~ 3 — 10~ 5 . After taking the upper bounds for the couplings, we checked their effects 
on the rare decay B® — > /J + /i~, the CP violation in D° — » h + h~ , and the isospin violation in B — > K^*' /j + /i~ . We 
found that the FC axigluon has virtually no effect on the decay B® — > since this process still occurs via loop 

diagrams. This result agrees with the last measurements of the branching ratio of this decay [2J, |25[. In the case of 
the isospin violation in B — !• A"'*)/i + /i - , FC axigluon effects seem to be at most the same order as the SM ones even 
though the axigluon contribution is at tree level. The most interesting effect of the FC axigluon is on CP violation 
in D° — > h + h~ decays. For this CP violating asymmetry, the upper bound on the axigluon contribution is at least 
10 times larger than the ~0.1% SM prediction. We conclude that adding small flavor violating effects to the light 
axigluon model might explain the CP violation in D a — > h + h~ and contribute to neutral meson mixings. 
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